A route to all frequency homogenization of periodic structures 
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We start from a one-dimensional periodic multilayered staclc in order to define a frequency power 
expansion of effective permittivity, permeability and bianisotropic parameters. It is shown from the 
first order that a simple dielectric multilayer can display a magnetoelectric coupling effect and, from 
the second order, that artificial magnetism can be obtained with arbitrary low contrast. However, 
this frequency power expansion is found to diverge at the first band gap edge. Thus, an alternative 
set of effective parameters, made of the propagation index and the surface impedance, is proposed. 
It is established that these effective parameters, as functions of the complex frequency, possess all the 
analytic properties required by the causality principle and passivity. Finally, we provide arguments 
to extend these results to the three-dimensional and frequency-dispersive case. 



There is currently a renewed interest in photonic crys- 
tals and metamaterials [Hll], i.e. periodic structures ex- 
hibiting new phenomena such as negative refraction |3JI3|. 
The former composites are known to possess ranges of 
frequencies (band gaps) for which no wave is allowed to 
propagate , while the latter composites exhibit a mag- 
netic response associated with local RLC-circuit type [S] 
or Mie [B], resonances. In this letter, we achieve such a 
magnetic activity with low-contrast dielectric layers. 

Our proposal is based upon a definition of all fre- 
quency homogenization (AFH) for which it is necessary 
to propose parameters satisfying causality principle [7l[8], 
i.e. analytic properties with respect to the complex fre- 
quency, and passivity. This necessitates to go beyond 
the Fresnel inversion [3] which, for each fixed frequency 
and wavevector, directly translates reflectivity, transmis- 
sion into effective permittivity, permeability and chiral- 
ity. The extension of classical homogenization theory 
[lOj to high frequencies [11) is of pressing importance 
for physicists working in the field of photonic crystals 
and metamaterials in order to understand extra-ordinary 
properties such as artificial optical activity [12] and mag- 
netism [51 [S] . Applied mathematicians show a keen inter- 
est in this topic [TTJ[T31[T3|, since periodic structures with 
small inductive and capacitive [5] elements structured 
at sub- wavelength length scales (typically A/10 to A/6) 
[316], can clearly be regarded as almost homogeneous. 

The starting point of this letter is the transfer matrix 
T(uj, k) associated with a layer which represents the unit 
cell of a periodic multilayered stack. The frequency uj 
and the wavevector k — (fci, fc2) are the conserved quan- 
tities of the system which is homogeneous with respect to 
time and space variables {xi,X2) [in this paper, we use 
an orthogonal set of coordinates {xi,X2,X3) such that 
the layers are stacked in the Xa-direction, see Fig. [4]. 
The transfer matrix T{z, k) is an analytic function of the 
complex frequency z (w is thus the real part of z) and the 
wavevector k. This analyticity property opens a route to 
the definition of effective homogeneous parameters for all 
frequency and wavevector spectrum. 

The matrix T(z, fc) is first used to build directly usual 



effective homogeneous parameters like the permittiv- 
ity ees{z,k), the permeability fj,ce{z,k) corresponding 
to artificial magnetism, and a bianisotropy coefficient 
^cs{z,k), hallmark of optical activity. The limit of this 
technique is then precisely defined: it is found that these 
effective parameters are appropriate only for frequencies 
ranging from the origin to the first band gap edge. 
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FIG. 1: Example of periodic multilayered stack of identical 
unit cells made of two dielectric homogeneous layers. The 
system is invariant under translations in the space directions 
x\ and X2 and the corresponding electric permittivity eix^) 
is a periodic function of variable x^. 

Next we turn to another set of effective parameters, 
the index ncs{z,k) and the impedance C,cG{z,k). Using 
general arguments based on existence of Bloch modes 
and the local density of states, it is shown that these pa- 
rameters are analytic functions in the upper half-plane 
of z, have a limit at infinite frequencies, and that imagi- 
nary part of ZTigg^z, k) and real part of Cos(z, fc) are both 
positive. As the central result of this paper, it follows 
that both effective index and impedance satisfy causal- 
ity and passivity requirements. This result leads us to 
the conclusion that a periodic multilayered stack can be 
replaced by a frequency- and spatially-dispersive homo- 
geneous effective medium for all the spectrum of frequen- 
cies w and wavevectors fc. The tool of choice for our one- 
dimensional model is the transfer matrix method, as it 
allows for analytical formulae, but we stress that ideas 
contained therein can be extended to frequency disper- 
sive and three-dimensional periodic structures. 

We first consider a periodic multilayer with a unit cell 
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made of two homogeneous layers L\ and £2 of thicknesses 
h\ and /12 ih = h\ + h2), see Fig. |4j At the oscillating 
frequency w, the electric and magnetic fields E and H 
are related to the electric and magnetic inductions D and 
B through the time-harmonic Maxwell's equations, 

-iujD{x) X H{x) , iujB{x) = V X E{x) , (1) 

and the phenomenological constitutive relations for non- 
magnetic isotropic dielectric media: 

D{x)^e^E{x), B{x) = fioH{x) , xeC^, (2) 

where fj,o is the vacuum permeability, and e„i the permit- 
tivity in the homogeneous layer Cm, rn — 1, 2. Note that, 
at this stage, the frequency dependence of the permittiv- 
ity is omitted. However, as it will be discussed later on, 
all the results of this letter remain valid when frequency 
dispersion is considered. In order to take advantage of 
invariance under translations in {xi,X2), a Fourier de- 
composition [from {xi,X2) to k = (fci,/c2)] is used to 
derive a first order differential equation from (|18|): 



dF 

- — (w, k, X3) = -iM,n{uj, k) F{uj, k, X3) , (3) 

where is a column vector containing the tangential 
components of the Fourier-transformed electromagnetic 
field {E,H), and Mm(cij,fe) is a matrix independent of 
x^ in each homogeneous layer For all vector x, 
is the component along the wavevector k — (fci, fc2) and 
x± its component along (—^2,^1)- Then, omitting the 
(w, fc)-dependence, one has for s-polarization 
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(4) 



Since the matrices are xa-independent, the solution 
of the equation ([s]) in each layer Cm is simply 



F{x3 + hm) = ex.p[-iMmhm] F{X3) . 



(5) 



The exponential above is well-defined as a power series 
of the matrix Mm, and defines the transfer matrix in the 
medium m through the distance hm ■ Since this power se- 
ries has infinite radius of convergence, the transfer matrix 
is analytic with respect to the three independent vari- 
ables id, fci and ^2 describing the whole complex plane. 
From now on, the complex frequency will be denoted by 
z = u! + ir], where lo remains the real frequency and 77 is 
the imaginary part. The transfer matrix associated with 
the unit cell, 

T{z, k) = exp[-iM2(2:, fc)/i2] exp[-iMi(z, k)hi] , (6) 

is also analytic function in the whole complex plane of 
variables z, fci and fc2 (for arbitrary permittivity profile, 
analyticity is proved using a Dyson expansion |15)). 



The analyticity property opens the possibility to ex- 
tract from the transfer matrix effective parameters which 
are analytic functions of the complex frequency z and 
valid over the whole frequency spectrum. First, the infi- 
nite radius of convergence of the power series expansion 
of T(z, k) suggests to introduce a notion of high order 
homogenization, which extends the usual homogeniza- 
tion (corresponding to the limit z — >■ 0) by expanding 
the effective permittivity and permeability as power se- 
ries of z. To carry out the asymptotic analysis, we use 
the Baker-Campbell-Hausdorff formula (BCH, extension 
of the Sophus Lie theorem, see [H]): 



exp[^] exp[S] = cxp[X] , 
X ^ A + B + lA,Bl + lA 



(7) 



where yl + _B is defined as the zeroth order approximation 
(classical homogenization), the commutator of A and B 
\A, B\ = {AB — BA)/2 is the first order approximation, 
lA — B, lA, BJl/3 is the second order approximation, and 
so forth. The BCH formula ([t]) shows that the transfer 
matrix ^ can be written as the one of a frequency- and 
spatially-dispersive homogeneous medium characterized 
hy X — — iMcff(z, k)h. The resulting matrix Mcs{z, k) 
corresponds to the constitutive equations 

D{k, X3) = Ecsiz, k) E{k, X3)+iKctfiz, k)JH{k, X3), 
B{k, X3) =fics{z, k)H{k, X3) + iJKcff{z, k)E{k, X3). 

(8) 

Here, matrix J represents the 90 degrees rotation around 
the xa-axis and, in the coordinate system (xy , x±, x^), the 
effective permittivity and permeability are 

Eoff = diag(e||,e^,e3), //off = diag(^|| , /i_L, Ma) , (9) 

while the bianisotropic parameter measuring the magne- 
toelectric coupling effect [121 Hi] is given by 



(10) 



i^eff = diag(X||,i^i,0). 



Although the system is isotropic in the plane {xi,X2), 
the spatial dispersion induced by the fc-dependence in- 
troduces a difference between the coefficients ey and e±, 
and /ij^, and K^^ and K± (TS]. At fc = 0, equalities 
£|| (z, 0) = e±{z, 0), and so forth. . . , are retrieved. 

After some elementary algebra, collecting terms up to 
the second order approximation in Q with A =^ — iA/2^2 
and B = —iMihi, we obtain the following homogenized 
coefficients [T^] for k = 0: 

£\\{z) = £1/1 +£2/2 + 52/i/2(ei-£2)(£i/i-£2/2)/(6eo), 
H\\{z) = Mo - z^/i/2Mo(£i-£2)(/i-/2)/(6eo), 
K\\{z) = i(ei-e2)/i/2/(2V£o7Mo), 

(11) 

with z = zh^/eoJiQ the normalized complex frequency 
and fm = hm/h. It is stressed that magnetoelectric 
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coupling comes from the odd order approximation in 
([7]), while artificial magnetism and high order correc- 
tions to permittivity emerge from even order approxima- 
tion. These results are fully consistent with descriptions 
in terms of spatial dispersion [T51 [5D] where, expanding 
the permittivity in power series of the wavevector, first 
order yields optical activity and second order magnetic 
response. This equivalence of these two descriptions (fre- 
quency and wavevector power series) is confirmed by con- 
sidering a unit cell with a center of symmetry, for exam- 
ple a stack of three homogeneous layers (permittivity 
and thickness hm, m — 1,2, 3) with £3 — ei and /13 = hi. 
Extending ([t]) to the case exp[A] exp[_B] exp[^] — exp[X], 
see [19j, it is found that K^s = 0, and thus it is retrieved 
that both magnetoelectric coupling and optical activity 
vanish in a medium with a center of symmetry [18j. 

Expansion in power series of frequency provides a new 
explanation for artificial magnetism and optical activity. 



Im(z) = r] 



"'^branch point 
-* — branch cut 



Rc(z) = uj 



FIG. 2: Domains of analyticity (gray regions) of effective pa- 
rameters Ecff, Mcff, Moff and ^off: the circle defines the disk 
where the power expansion converges. 



Keif — 0: omitting the (2, A;)-dependence, we have 



Tr 



off ■ 



cos[zncBh] — ^(Ccff/^:) sin[znoff/i] 

-i(z/Ccff) sin[zncff/i] cos[zncsh] 



, (13) 



where z^nlg = z^EoffMcff - k'^ and Ccff = Moff/"-cff- Com- 



Analytic expressions (11) of effective parameters can be parison of (12) and (13) gives definitions of propagation 



used to analyze artificial properties (expressions for nor- 
mal incidence up to order 6 are provided in |19j). In 



index (along xa-axis) and surface impedance Ccff: 



particular, we show from (11) that: artificial magnetism. 



previously proposed with high contrast |B1 E] , can be 
obtained with arbitrary low contrast; and magnetoelec- 
tric coupling, previously achieved in fi-composites |12j . 
can be present in simple one-dimensional mulitlayers. 

Nevertheless, this frequency expansion of effective pa- 
rameters cannot be used for frequencies higher than the 
frequency wi at the first band gap edge ^9j. To show 
this limitation, we consider for the sake of simplicity a 
3 layer unit cell with a center of symmetry and purely 
real dielectric constants ei = S3 and 62- effective pa- 
rameters EcS and /ieff derived from ([t]) are then purely 
real and K^s = 0. At the band gap edge uji, either EoS 
or /Lteff must vanish to allow for sign shifting, while the 
other parameter must take an infinite value: the power 
series expansion of X = —iMcgh in ([?]) diverges at uJi. 
Indeed, when calculating X — log{exp[y4] exp[i3]}, it ap- 
pears that the function log is not analytic when its ar- 
gument "vanishes" , which introduces a branch point at 
uJi: this branch point implies that the radius of conver- 
gence of the power expansion is bounded by loi (see Fig. 
[5]). However, it is possible to choose the branch cut of 
the complex logarithm from the branch point wi in the 
lower half complex plane (Fig. [5]). This choice makes 
it possible to define effective parameters which are an- 
alytic functions in the upper half plane of the complex 
frequency z, as required by causality principle |18j . 

We consider a unit cell with a center of symmetry to 
keep things simple (no magnetoelectric coupling). The 
general expression of the transfer matrix is then [2T] 



cos[znoe{z, k)h] = a{z, k) , Ceff(2, fe) = 



lb{z,k) 
d{z, k) 



(14) 

Here, notice that zucsiz, k) is just the Bloch wavevector. 
The main result of this letter is the following. In the 
upper half plane of z, i.e. for Im{z) > 0; i) imaginary 
part of zneff{z, k) and real part of Ceff{z, k) are positive; 
ii) neff(z,k) and C,eff{z,k) are analytic functions of z; 
Hi) n^fflz, k) and Ceff{z, k) have limits n^o and Coo when 
\z\ — > 00, where n^o = {^EfJ-o ) is the mean index. 

To prove these claims, we first use the theorem stating 
that no Bloch mode exists for z in the upper half plane 
[22j . As a consequence, the function zncs{z, k) cannot be 
purely real and its imaginary part Im(znoff) cannot van- 
ish. This proves assertion i) for Uos. Next, it is stressed 
that the coefficient a(z, fe) is z-analytic in all the complex 



plane of z, and that the definition ( 14 ) can be written 



exp[izncs(z, k)h] — a{z, k) + i\/\ — a'^{z, k) . (15) 

Since Im(znoff) cannot vanish if Im(z) > 0, we have 
cos[2;noff (z, k)h] = a{z, k) ^ ±1, and thus the square root 



in ( 15 ) is z-analytic in the upper half plane. The func- 



tion on the left hand side of ( 15 ) is then analytic and. 



in addition, cannot vanish. The complex logarithm can 



be applied to ( 15 ) without alteration of the analyticity 



property: this proves ii) for n^s. Next, combining the two 



equations a —bd — 1 ( 12 ) and a ^ ±1 for Im(z) > 0, it is 



T(z,fe) 



a{z,k) b{z,k) 
d{z,k) a{z,k) 



bd^l. 



(12) 



found that none of the two analytic functions b{z, k) and 
d{z,k) vanishes, and thus the ratio b{z , k) / d{z , k) ^ 
is analytic. The square root in (14) preserves the ana- 



This matrix is compared with Tcsiz, fe), the transfer ma- 



trix corresponding to the constitutive equations (25 ) with 



lyticity property, which proves assertion ii) for Ceff- The 
proof of i) for Ccff is based on the local density of states. 
The Green's function of the multilayer is calculated [53] 
with a point source located in the plane ^3 = : the 
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value of the electric field in the same plane is found to be 
*7''Ccff- Since the imaginary part of the Green's function is 
positive (it corresponds to the local density of states), it 
follows that Ccff has its real part Re(Ceff) positive. Finally, 
the proof of iii) is given in [TH] . 

The main result tells us that artificial frequency dis- 
persion, i.e. z-dependence (of effective parameters) gen- 
erated by periodic spatial modulation, has the same prop- 
erties as the natural frequency dispersion in usual media. 
Part i) ensures passivity, while parts ii) and iii) imply 
that the effective parameters fulfills causality principle 
[THj . From ii) and iii) , the Cauchy integral formula can be 
applied to the function ncg{z, k) — rioo to obtain Kramers 
and Kronig relations and their generalization ^ : 

... 1 / , J/lm[neff(i^,fc)-noo] 

ncs(z,k) = noo / aiy ^ ^ . (16) 

Jb z — v 

An illustrative example in Fig. |6]confirms that this gener- 
alization of the Kramers and Kronig relations is satisfied 
by nefi(-2;, fc), since the solid curve and plus markers of 
Re(neff) fit each other. Also, Fig. 3 confirms that, at the 
infinite frequency limit, the effective index tends to the 
mean refractive index (y/e), as in usual media for very 
high frequencies of neutrons [23]. Note that, contrary to 
frequency independent dieletric constants £i = £3 = 2eo 
and £2 = 12eo of the multilayered stack, effective param- 
eter noff(z, k) agrees with causality principle. 

The main result remains valid when (natural) fre- 
quency dispersion is considered in the dielectric permit- 
tivity [e.g. dielectric constants Emiz)]. Indeed, permit- 
tivity is analytic in the upper half plane of z and theo- 
rems on existence of Bloch modes and imaginary part of 
Green's function can be applied in the most general cases 
[TIES]- Note that, when natural frequency dispersion is 
considered, the limits n^o and take the values of index 
and impedance in vaccum, i.e. y/eoHo and -y^/io/^o- 

Finally, the possibility to extend our results to fre- 
quency dispersive three-dimensional periodic structures 
is discussed. Using the auxihary field formalism [7] l22]. 
Maxwell's equations can be written as the unitary time 
evolution equation [dF/dt]{t) = —iKF{t), where K is self- 
adjoint and time-independent (see also the frame of ex- 
tension of dissipative operators [25]). Consequently the 
resolvent [z — K]~^ is analytic if Im(z) > 0, which pre- 
vents the existence of Bloch modes in the upper half plane 
of z. And imaginary part of Green's function is always 
positive for lm(z) > since the operator 



1 



1 



K 



K 



lm(z) 



1 



1 



Kz-K 



> (17) 



is positive. Assuming that appropriate effective parame- 
ters can be defined (see promising notions of impedance 
in [5BH2H]), these two properties on Bloch modes and 
imaginary part of Green's function can be used to prove 
causality principle and passivity. It is however stressed 
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FIG. 3: Real and imaginary parts (solid lines) of effective 
refractive index ncfi{z) for z — uj + 0.001 x i, deduced from 
(151; Kramers and Kronig formula unveils Re(ncff) in plus 
markers; ei = £3 = 2eo, £2 = 12£o, /i = /a = 0.4, /2 = 0.2, 
\/{£) = Vei/i + £2/2 + £3/3 = 2eo, {V^) = V£ifi + V^f2 + 
V^/3 ~ 1.82eo. 



that such a generalization remains a challenging task: 
particularly, special attention should be paid to sit- 
uations where single mode Bloch approximation does 
not apply [28,, and to analytic continuations of Bloch 
wavevector and impedance for complex z. 

In conclusion, a periodic multilayered stack can be 
modelled at low frequencies by a homogeneous medium 
characterized by effective permittivity, permeability and 
magnetoelectric coupling parameter. But these param- 
eters, defined as power expansions, are not approriate 
for frequencies beyond the first band gap edge. Con- 
sidering instead effective propagation index and surface 
impedance, it is shown that artificial frequency dispersion 
has the same properties as natural dispersion in terms of 
passivity and causality: remarkably it follows that a pe- 
riodic arrangement of frequency independent (and thus 
non-causal) dielectric materials makes artificial causality. 
These results, based on general properties of existence of 
Bloch modes and sign of imaginary part of Green's func- 
tion, open a route for generalization of AFH to frequency 
dispersive and three-dimensional case. 
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SUPPLEMENTAL MATERIAL 



Effective parameters for a one-dimensional periodic multilayer 



In this supplementary material we explain in more details our homogenization algorithm for the periodic multilayer 
as shown in Fig. |4j The effective medium is depicted in the right panel, with effective permittivity, permeability and 
bianisotropic parameters being rank-2 tensors. Let us start from the Maxwell's equations 
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(a) (b) 
FIG. 4: Schematic diagram of the periodic multilayer and corresponding effective medium. 



- iujD{x) = V X H{x) , iujB{x) = V X E{x) , (18) 

and the constitutive relations for non-magnetic isotropic dielectrics 

D{x)^ernE{x), B{x)=^loH{x), X e C^, (19) 

where /xq is the vacuum permeability and e,„ the permittivity in homogeneous layer Cm- Here we introduce the 
Fourier decomposition for both the electric and magnetic field in the form 

E(fci, fc2,X3)= / 'E{xi,X2,X3)exp[-i{kiXi + k2X2)]dxidx2, 



H(fci, ^2, xa)^ — / li{xi,X2,X3)exp\ — i{kixi -\- k2X2)]dxidx2 ■ 



(20) 
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and a 4-coniponents column vector [1^: 



Applying the decomposition (20) to (18 1 and (19), we derive an ordinary differential equation involving 4 x 4-matrices 

(w, k, xs) = -iMm{i^, k) F{oj, k, xs) . 



dF 



dx-. 



(21) 



Here, the field-components and H3 have been eliminated and the resulting column vector F contains the electric 
and magnetic field-components which are tangential to the interface of the multilayered stack, i.e. the components 
along xi and X2- 

In order to simplify the derivation, a coordinates rotation of the (xi, a;2)-plane is used: The new coordinates are 
denoted by x\\ and x± where the component xy is along the wavevector k — (fci,fc2) and xj_ is along (— /c2,fci)- In 
other words, the new coordinates are related to the previous ones through the following transformation 



x± 



1 



2 I l2 
1^2 



ki k2 




Xi 


~k2 ki 




X2 _ 



Thus, the column vector F for a s-polarization incidence is 



Correspondingly, the matrix is 



ujH\\ 





,2^ 7.2 , 



Mo 




LO^£,n - k 

For the effective medium in the right panel of Fig. |4j the constitutive relations turn out to be 

Dik, x^) = ecff(w, fc) E{k, xz)+iKcf[{(jj, k)JH{k, X3), 
B{k, X3) = /icff(a;, k)H{k, X3)+iJKcs{uj, k)E{k, X3), 
where egg and jieS are the effective permittivity and permeability tensors, 



(22) 



(23) 



(24) 



(25) 





" ^11 





" 









" 


Eeff — 






















. 





£3 . 




. 





Ai3 . 



(26) 



Matrix J represents the 90 degrees rotation around the xa-axis, and KcS is the bianisotropic parameter measuring 
the magnetoelectric coupling effect: 





"0 


-1 


0" 




' 





0" 


J = 


1 








, Kcft ~ 












. 





1 _ 














(27) 



Similarly, we can obtain the matrix M^s by adopting the same calculation process of M„j in a periodic multilayer 

—iujKu /ill 



UJ £j_ — fc //i3 iw/i'ii 

The exponential function of matrix Mm defines the transfer matrix for the multilayer 

T{uj, k) = cxp[— iAf„j/i„j] 

while for the effective medium in the whole complex plane, its transfer matrix is 

Tcff(i:j, k) ^ cxp[-iAfoff/i] ■ 

They should satisfy the following relation: 

exp[— jMcff/i] = exp[— 1^2/12] exp[—iAfi/ii] . 



(28) 

(29) 
(30) 
(31) 
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In order to solve this equation, we need to introduce the Baker-Campbell-HausdorfF (BCH) formula '5], which defines 
an approximation for the product of two exponential functions with matrix as arguments 

exp[yl] exp[B] = exp{A + B + lA, Bj + IjB - A, IB, A]] + • • • ) (32) 

here we denote A + B the zeroth order approximation, the commutator JA, BJ = {AB — BA) /2 the first order 
approximation, and so forth. Hence taking A — —iM2hi, B — —iM\h\, we obtain 

-zMeff/i = ~t{M2h2 + Nhhi) + I-^M2/^2, -iMi/iil + ^hiAf2/i2 + iMihi, [-zAf2/i2, -iATi/ii]] + • • • (33) 



Assuming a normal incidence whereby fe = 0, and taking the zeroth order approximation in (33), we have 

Afeff « Afi/i + Af2/2 (34) 

with /i = hi/h, ji = hijh, which is the result the classical homogenization gives [3]: 

£_L = £i/i + £2/2, Mil = Mo, K\\ = . (35) 

One should note that classical homogenization does not capture any artificial magnetism or bianisotropy effect. 
Next, taking up to the first order approximation in (33) we obtain 

A/eff « Afi/i + A/2/2 - ilM2/2, ATi/il , (36) 
and the effective parameters turn out to be 

ej_(z) = £1/1 +£2/2, M||(i)=Mo, i^||(z) = z(£i - e2)/i/2/2eo , (37) 

which produces a magnetoelectric coupling parameter K depending on the frequency: here, z = zh^/eoiiQ is the 
normalized complex frequency (the complex frequency is z = lo + irj, its real part being the real frequency lo). 
The second order approximation in ( 33 1 yields 

Afeff « (Af2/2 + Afi/i) - zlA/2/2, Mifij + ^pf2/2 ^ A/i/i, I-^A^2/2, -*A/i/ill. (38) 

Based on this, we can finally obtain the analytic expressions for effective parameters in equation (11) of the letter, 
where the artificial magnetism and high order corrections to permittivity appear. It should be noted that, under a 
normal incidence, the s-polarization and p-polarization waves are the same under the isotropic system in the plane 
{xi,X2), i.e. e||(z) = s±{z), ^i\\{z) = ^J,±{z) and K^iiz) = Kj_{z). 

Based on the BCH formula, we can push our algorithm to any order approximation. Here we consider up to the 
sixth order approximation where the effective parameters are 

eu{z) = ej_iz) = £1/1 +£2/2 + -^/i/2(£i - £2)(£i/i - £2/2) 

II ^ 

- ^/l/2(£l - ^2) [(£2/1 - £l/f )(£l/l + £2/2) + 3/i/2(£l/2 - £l/l)] 

+ ^/i/2(£i - e2)(4£l/f + 18£1AV2 + 27£l/f/| + 6e?e2/iV2 + Seje^fHi + Se^fHI 

- 3£i£l/iV2 - Sii4f!f2 - 6£i£^/i/2' - 27£l/2/| - 18£l/i/| - 4£1/|) , 
Mll(^) = M-l(^) = Mo - yMo/i/2(£i - £2)(/i - /2) 

+ ^M0/l/2(£l - £2) [(£2/1 - £l/f ) + 3/i/2(£l/2 - £2/1)] 



^Mo/i/2(£i - £2)(4£f/f + 6ef /f /2 - Sefftf^ - 21eif(f^ + 1^162^2 + ^Wtfi 
8£i£2/f /I - 18£i£2/i/2" + 27e2/3/2 + ulflfl - 6£l/i/| - Aelfl) , 



X|| (z) = KAz) = 2 (£1 - £2)/i/2 + el)flfl 



^/l /2 (£1 - £2)(£?/r + 3e1/i/2 + £l£2A' + 2eie2/l/2 + £l£2/2' + 3£l/l/2 + £1/|) , (39) 
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where ii = Ei/eo is the relative permittivity (eq is the vacuum permittivity). It is clear from this expression that all 
effective parameters depend crucially upon frequency and, permittivity and filling fraction contrasts. 




FIG. 5; Effective parameters versus normalized real frequency 



The curves of these effective parameters versus normalized frequency are plotted in Fig. [5] where we have assumed 
the two dielectrics in the unit cell of the periodic multilayer are Glass and Silicon respectively, with 

£1 = 2, £2 = 12, /i=0.8, /2-O.2. (40) 

It can be observed that all the effective parameters increase along with the frequency, the effective permittivity takes 
the values larger than 1 (artificial magnetism), while K describing the magnetoelectric coupling is always non-zero. 
One should note that all correcting terms in £||(z)[=: £_l(z)] and ^\\{z)[— /ij^(z)] are positive. 



BCH formula for the 3 layers' case 



For a periodic multilayer with a unit cell consisting 3 dielectric layers, the relation between the transfer matrices 
of multilayers and its effective medium will be 

exp[X] = exp[A] exp[B] exp[A] , (41) 

which is a product of three exponential functions with matrix argument. In order to derive the expressions for X , we 
can take an iteration of BCH formula of (32) using 

exp[C] = exp[A] exp[i?] 

with 

C = A + B+\lA,B\ + ^lA,lA,B\\-^ {B, {A, B}}--- 



(42) 



(43) 



Then equation ( 41 ) becomes 



exp[X] = cxp[C] exp[^] , 



(44) 



and hence 



X = A + B + A + ]^IA,B\ + ]^IA + B,A\ + ^IAM.B\\-^ IB, {A, B\\ 
+ J [[A, BlAl + ^^lA + BM + B, All - ^ lA + S, All + . . . 



(45) 
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Similarly, we take up to the fourth order approximation under a normal incident in s-polarization, the effective 
parameters are 



e||(z) = ej_{z) = 2£i/i +£2/2 - y/i/2(£i - e2)(ei/i +£2/2) 



45 



/i/2(£i - i2)i7eift + 3eiftf2 + ll£i£2///2 + 2£i£2/i/2' + iifi + Qeihfi) , 



/i||(z) = ^j_{z) = + yAio/i/2(ei - £2)ifl + /2) 



45 



Mo/i/2(ei - e2)(7£i/f + ll£i/f /2 + 6£i/i/| + '3e2f!f2 + £2/! + 2£2/i/|) 



K^^{z) = Kj_{z) = . (46) 

The results indicate that the magnetodielectric coupling vanishes in the three layers' case, which can be explained 
by the fact that: the odd order approximations in (45 1 are all equal to zero due to the anti-commutation in the 
commutators ([[A, B] = -lB,AJ). 



Divergence of the power expansion at the band gap edge 



In 



this section, we show that expressions of the effective parameters introduced by the power expansion in previous 
sections are no longer valid at the first band gap edge. 

pand the transfer matrix in (30) with Taylor series: 



actions are no longer valid at the hrst band gi 
First, we expand the transfer matrix in (30) 



T,ff = exp(~^A^eff/^) = ^(-^)2^M_ + ^(_,)2p+ii^ 



^^2p+1^2p+l 



p=0 



According to (28) where the real frequency uj is replaced by the complex one z, we 



(2p+l)! 
have 



— IZ 
2 



zi^ll/co 



Mil 



- / IJ.3 izK\\ /cq 



V 4 " y 



1 
1 



" '''cff I 



with 



Ai3 



Let us plug (48) into (47) and, considering Taylor series of the sin — cos functions 

sin(A) = y / ^ , , cos(v4) = y ^-^^2" 

n=0 V ' ; 



•1=0 



(2n) 



we obtain 



Teff = cos{kcsh) - i^^^^ sin(fcoff/i) = 

Kcff 



COs(fcefl/l) 

-i(z2ej_ 



zifjl sin(A;eff/i) sin(fccff/i) 



"7 

Co Kcff 

fc^ sin(fcoffM ,^ 
— ) cos(fccfi /i) - 

Koff 



fceff 

zi^il sin(A:cff 



L f^d '"ott Cg fccff 

Let us now write the characteristic matrix |3] of the dielectrics of the periodic multilayer as follows 



COS (kihi) sin{kihi) 

z ^ 

— z — sm{kihi) cos{kihi) 

Ci 



ki — Z Cifii A? , — 



(47) 



(48) 



(49) 



(50) 



(51) 



(52) 



The transfer matrix of the unit cell consisting of three dielectric layers (ei = £3, hi = /13) is 

T = T1T2T1 . 



(53) 



10 



Since the dispersion law is defined by the trace of transfer matrix T, we have 



tr(r)/2 



COs(fcoff /l) 

COS (2fci/ii) COS (^2/12) — 0(7^ + 7^) ^™ (2fci/ii) sin (fc2/i2 
2 Ci C2 



effective medium, 

muhilayer stack with three layers. 



(54) 



The dispersion law of the effective medium is depicted in dashed line [here we have taken up to the 20th order 
approximation in (45), and substitute the analytic expressions of the effective parameters into (49) and then to (54)], 



as well as that of the multilayers in solid line in Fig. |6] assuming a normal incidence within s-polarization. One 
can see that the dashed line diverges when frequency tends to the lower edge of the first stop band, Sdenoted by 
LJi. In other words, the expressions of these effective parameters are not valid any more further the first band gap in 
multilayers. 




FIG. 6: Dispersion laws of multilayers (solid line) and effective medium (dashed line in 20th order approximation). 



Proof of iii) of the main result: Effective index and surface impedance at infinite frequency 

As it has been shown in the previous section, the effective permittivity, permeability and magnetoelectric parameters 
are no longer valid when uj > uji. To overcome this pitfall, we consider the set of effective parameters made of 
propagation index and surface impedance Ccff- 

We consider a periodic multilayer with a unit cell consisting of three layers, and comparing T = T1T2T1 with the 
transfer matrix of the effective medium in (51), we find 



and 



1 C2 Ci 

a = cos [zricfih] = cos (2fcift-i) cos (^2/12) — o + 7^) ^™ sin (fc2^2) 

? ^1 > ^"^ 

= cos (2fcift,i + fc2/),2) - ;:,(7^ + 7^ - 2) sin {2kihi) sin (fc2/i2) , 
2 Ci C2 

C 2 2 

zh — —i2(i sin (kihi) cos (fci/ii) cos (A;2^2) + i-r~ sin (fci/ii) sin (fc2ft.2) — iC2 cos (/3i) sin (^2^2) , 

2 1 

d/z = —i— sin (fci/ii) cos {kihi) cos (^2/12) + sin {kihif' sin (^2/12) ~ cos (fcift-i)^ sin (^2/12) 
Ci Ci C2 

When frequency z = uj + irj — |z|e*'* tends to infinity by \z\ — )■ 00, Euler's formulae 



simplify into 



(55) 



(56) 



sinz= , cos 2;= , (57) 

2i 2 ^ ^ 



e e 

sinz~ , cosz~ , (58) 

2t ' 2 ^ ^ 
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since e*^ — >■ 0. Equation (55) becomes 



exp {—izucsh) ~ exp {—i2kihi — ik2h2) + Aexp (— i2fci/ii — ifc2/i2) 
with A = (C2/C1 + C1/C2 " 2)/4 a constant. Furthermore, 

izuesh ~ i{2kihi + A;2/i2) — In (1 + A) 

when |z| — >■ 00, — >■ 2:71^, then 

rioff (2zni/ii + zn2h2)/{zh) H ln[l + A] ^ (n) . 

zn 00 



(59) 
(60) 

(61) 



This proves that the Umit of the refractive index at infinite frequency is equal to the mean of indices of a unit cell, as 
stated in the letter. 

The surface impedance is defined by 



^2 _ 2 

''Off — ^ 1 



2 2 
-i2(^i sin (fci/ii) cos {kihi) cos (fc2/i2) + i-r- sin (kihi) sin (^2^-2) ~ iC2 cos (/3i) sin (fc2/i2) 

C2 

d 2 C2 2 f 2 

—I— sin [kihi] cos (kihi) cos (^2/12) + *72 ^™ (^1^*1) sin (A:2ft-2) — cos (kihi) sin (^2/12) 
Ci Ci C2 

sin (fc2/i2) [Ci sin^ (kihi) — Q cos^ (fci/ii)] — (1C2 sin (2A:i/ii) cos (^2^2) 



sin (fc2/i2) [Cl sin^ (kihi) — cos^ (fci/ii)] — C1C2 sin {2kihi) cos (fc2/i2) 



(62) 



Plugging (58) into (62), 



Coo = lim Ccff = lim z-y/ - 



lim < 

\z\^oc 





'a 


_g2ifci/ii 


-Cl 




-C1C2 


g2ifci/ii gi 




2i 


4 


4 


2i 


2 




[a 


^2iki hi 






-C1C2 




k2h2 


2i 


4 


4 


2i 


2 



1/2 



(63) 



the surface impedance is shown to converge at infinity to the impedance of layer 1. 

More generally, for a unit cell consisting of m layers with permittivity e,„ and thickness /i„i, the transfer matrix is 



(64) 



When \z\ — > 00, (52 1 becomes 



exp [—ikihi) 



1 -K[ 
1 



with C,[ ~ C,i/ z [i ~ 1, . . . m) . According to (64), we have 



Til ^ exp [-iz{kihi + fc2/l2 h fcmft.m)](l + ^) 

where A is an algebraic function of ratios CiKj — CUCj {hJ ^ ^7 ■ ■ ■^)- Finally, we obtain 

exp (—izricsh) ^ exp [~iz{kihi + A;2^2 • • • + fcm/Jm)](l + A) , 

and 

n-csiz) ^ (kihi + k2h2 h k.mhm)/{zh) + — ln[l + A] ^ (n) . 

zn 00 



(65) 

(66) 
(67) 
(68) 



The same limit for the effective parameters nce{z) and Ccs{z) can be obtain when Re(z) = uj — > 00. In this case, 
adding an arbitrary small imaginary part to dielectric constants Em, the above derivation applies mutatis mutandis. 
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Kramers and Kronig formula for the effective surface impedance 

From the definition of effective surface impedance given by equation (14) in the letter, we know that the function 



F{z,k) = i[Ccs{z,k) - Co 



(69) 



is analytic in the upper half-plane, and tends to when \z\ — > oo. Applying the Cauchy integral formula to function 
F{z), a relationship equivalent to Kramers and Kronig formula can be obtained: 



[i(Ceff(2,fe) - Coo)] = 



1 



dv 



Im [i(Coff(j^, fc) - Co 



(70) 



The real and imaginary parts of F{z) are shown in Fig. [t] in solid lines, while the line denoted by plus markers is 
obtained from ( 70 ) , while a small dissipation has been introduced in the permittivities of the dielectrics to enforce 
the convergence of the effective surface impedance. The consistency between those two curves of the real part of F{z) 
confirms that Ceff must satisfy the Kramers and Kronig formula. 





FIG. 7: Real and imaginary parts (solid lines) of effective surface impedance i{(^ca — Coo) for z — cu + 0.1 * i, deduced from 
(70l; while Cauchy integral formula unveils Re [^'(z)] in plus markers; here ei = £3 = (2 + 0.1 * i)eo, £2 = (12 4- 0.1 * i)£Q, 
h^h^ 0.4, /2 = 0.2. 
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